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Abstract. We introduce a class of one-dimensional lattice models in which a 
quantity, that may be thought of as an energy, is either transported from one 
site to a neighbouring one, or locally dissipated. Transport is controlled by 
a continuous bias parameter q, which allows us to study symmetric as well as 
asymmetric cases. We derive sufficient conditions for the factorization of the N- 
body stationary distribution and give an explicit solution for the latter, before 
briefly discussing physically relevant situations. 



PACS numbers: 05.40.-a, 02.50.Fy, 47.27.Eq 

Introduction 

Systems driven in a non-equilibrium steady state by an external forcing that generates 
an internal flux are very frequent in nature. Among these systems, two broad classes 
may be distinguished according to their conservation properties. The first class 
corresponds to systems in which the flux of a conserved quantity (e.g., particles 
or mass) takes place. Such systems have raised considerable interest recently, and 
different paradigmatic models for which the full (./V-body) stationary probability 
distribution can be worked out exactly have emerged from these studies: the 
Asymmetric Simple Exclusion Process (ASEP) ^ 12 El ED H] , the Zero Range Process 
(ZRP) PEHH, the Asymmetric Random Average Process (ARAP) El EH 03] , as 
well as more general mass transport models 1121 . 

On the other side, a second class corresponds to situations where the quantity 
moving through the system is locally non conserved in the bulk. This may happen 
for instance when the system can exchange particles with a reservoir, or when one 
considers an energy flux in a system where the 'microscopic' dynamics -at the chosen 
level of description- is already dissipative, like in turbulent flows or in shaken granular 
materials for instance. Within this class of systems, attention has been mainly devoted 
to ASEP or ZRP models where particles can be added or removed within the bulk, at 
rates that may differ from the boundary ones |7l 1141 IT51 1161 1171 118) . Yet, there seems 
to be very few known solvable models, for which the probability of any microscopic 
configuration can be computed exactly, where a continuous quantity is injected at the 
boundaries and (partially) dissipated in the bulk, a situation of broad physical interest 

cat- 

In this note, we consider a generalisation of the one-dimensional cascade model 
introduced in |20| . where a local quantity, injected at the boundaries, is either 

X Note that a 'reverse' situation, where a quantity is constantly added in the bulk, has been studied 
in the context of force fluctuations in beads packs, the vertical axis playing the role of time 1101 . 
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Figure 1. Schematic view of the model, with the three different mechanisms: 
energy injection at boundary sites with rates J^(/ii) and Jfl(Mjv), transport from 
site i to site i + 1 and i — 1 with rates q<f>i(fij\pi) and (1 — q)(pi{pi\pi) respectively, 
and dissipation at site j with rate Aj(/j,j\pj). On sites i = 1 and i = N, energy 
may also be transferred to the reservoirs —see text. 



transported or dissipated in the bulk. A bias parameter q allows us to consider partially 
asymmetric transport. We give an explicit solution for the full ./V-body stationary 
distribution in cases where it factorizes, and briefly discuss a physical application to 
a schematic turbulence modeling. 

Definition of the model 

We consider a one-dimensional lattice with sites labelled from i = 1 to N, and 
introduce a local quantity pi, that may be thought of as an energy. The dynamics is 
defined through the following local stochastic rules. An amount of energy between p 
and ji + dp can be injected from a reservoir on the left boundary i = 1 with a rate 
(probability per unit time) Jj 1 {p)dp, or injected from another reservoir on the right 
boundary i = N with a rate Jn(fi)dfi. Transport proceeds by removing an amount of 
energy /i on site i with a rate 4>i(n\pi)d^, depending on the energy pi present on site i. 
The energy p is then added either to site i + 1 (to the right reservoir if i = N) with 
probability q, or to site i — 1 (to the left reservoir if i = 1) with the complementary 
probability p = 1 — q. Finally, dissipation is implemented by withdrawing an energy p 
from site i with a rate Ai(p\pi)dp, without adding it to another site. A simple sketch 
of the model is shown on Fig. ^ 

The statistical state of the model is described by the TV-body probability 
distribution P(pi, pN,t)- Given the above stochastic rules, the evolution of this 
distribution is governed by the following master equation: 

dP v— \ f Pj+1 

— / dpcjjjiplpj + p)P(...,pj + p,p j+ i- p,...,t) 

Ct j=1 Jo 



+P2_, / dp(j)j(p\pj + p)P(...,pj-! - p,pj + p,...,t) 
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where the dots stand for all the variables P i that are not modified by p. 
Factorized steady state 

We wish to study the stationary distribution P s t{pi, Pn) in cases where it is 
factorized, and thus make the following ansatz for P st ( P i, ..., P n)- 

1 N 

P s t(pi,...,PN)=zl[fi(p i ) (2) 

i=l 

where Z is a normalisation factor. Note that due to the presence of dissipation in the 
bulk, one expects that the one-site distribution fi(pi) is different from one site to the 
other, even when the stochastic rules are site-independent, i.e., <pi(p\pi) — (f>(p\pi) 
and Ai(p\ Pi ) — A(p\ Pi ). Introducing the above ansatz for P st (pi, p N ) in the 
stationary master equation, it is easy to check that a specific form for the transport 
and dissipation functions (f>i(p\pi) and Ai(p\pi) allows for a decoupling of the different 
degrees of freedom, in the sense that the stationary master equation can then be 
written as a sum of terms involving at most one variable pi. This specific form reads: 

H»\p i ) = Uti fi{P f i ~t b \ A i (p\p i ) = A i (p) fi ^ i - lx) (3) 
TAPi) TiKPi) 

Under these assumptions, the stationary master equation can be rewritten as: 

Y r^fe(M)+A j ( M )] /j( / ? ;~ M) + r dp[j L {p)+j R {p)\ 

Jo J]\Pj) Jo 
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and the decoupling property mentioned above appears explicitly. As this equation is 
a sum of functions of independent variables, each of these functions -i.e., the sum of 
terms depending on a given p^— has to be equal to a constant. To determine the value 
of these constants, one can first send all the p\s to zero, in which case all the integrals 
over the interval [0,Pi] are expected to vanish §. As a result, the terms that do not 
depend on pi sum up to zero: 

J o dp, Lfo(ri+pfa(p)+J2*i(»)- Jl(p)- Mv)j =0 (5) 

In a similar way, sending to zero all the pis but one leads to the conclusion that all 
terms depending on a given pj also sum up to zero. For 2 < j < N — 1, one has: 

P3 dp [qfc-iQt) + pfa+i (p) - fa (p) - Aj (p)} !&LzA = (6) 
o JjKPj) 

At the boundary sites j = 1 and j = N, one finds: 

P1 dp [J L (p) + pfa (p) - fa (p) - A x (p)} fl{pi - M) = (7) 
o h(Pi) 

Pl dp [J R (p) + q^N-^p) - 4> N (p) - A^(/i)] /jv(P f ~ M) = (8) 

o Jn(Pn) 
A sufficient condition to solve this set of equations is to assume that each of the 
integrands is equal to zero, leading to: 

JV 

q4> N (p) +PMP) + A i(<") _ - Jr ^) = ( 9 ) 
3=1 

qfa„ 1 (p)+pfa+ 1 (p)-fa(p)-A 3 (p) = (2<j<N-l) (10) 

J L {p)+p4> 2 {p)-4> 1 {p)-h 1 {p) = Q (11) 

J R (p) + q<t> N -i(p) - 4>n(p) - Ajv(m) = (12) 

Thus, as a consequence of the specific form @ chosen for <fij(p\p) and Aj(p\p), it turns 
out that the function fj{p) has disappeared from the equations -although this result 
may be surprising at first sight, it should be noticed that the mass transport model 
introduced in \T2\ has the same property. As a result, any set of functions fj{p) may 
be a solution of the model, provided that the dynamical rules are suitably chosen. 
Eqs. ©, (fTTTf) . (fTTJl and (fT2"|l should thus be considered as compatibility conditions 
allowing for the existence of a factorized distribution. Let us also note that these four 
equations are not independent: it can be checked easily that, for instance, Eq. @ can 
be obtained from the sum of Eqs. (|10fl -summed over j—, 1|11|) and (|12fl . To solve these 
equations, we first define some extra functions (f>o(p) and 0jv+i(/x) on the boundaries 
as -assuming p, q ^ 0: 

J L (p) = q4> (p), J R {p) =p4> N+1 (p) (13) 

which allows one to rewrite Eqs. (|ll|l and (|12|l on the same form as Eq. I|10J) with 
j = 1 and j — N respectively. Introducing an auxiliary function Xj (/■*) through: 

X] (p)=fa+i(p)-h(p), 0<j<N (14) 

§ Considering an integral of the form I(p) = d/j,il>(fj,)f(p — (J.)/f(p) with V(m) ~ fJ. a ~ 1 and 
f(p) ~ pfi" 1 for fj, — * (a, P > to ensure the convergence of the integral), it can be shown easily 
that I(p) ~ p a -» when p -> 0. 
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Eqs. QlOfl. |(TJ| and (|12[l can be rewritten in a concise form as: 

PXj(p) - <%-i(m) = A j(m)> 1 < 3 < N (15) 
Setting r = p/q, this last equation is easily integrated out as: 



N-l 

xM=K(p)r N -i --£V-'A i+1 0<j<AT (16) 



where K(p) is (up to now) an arbitrary function of /i, and with the convention that 
the sum is zero if j = N. Using Eq. (|14fl and assuming q ^ 1/2, one can now solve for 
<f>j(n), yielding for < j < N + 1 -the sum is zero if j > N: 

1 N 

4>j (m) = ^(m)^-^ + — - J2 (i - ^')A 2 + <?(m) ( 17 ) 
9 t=j+i 

with a new (arbitrary) function C(/i), and where K(p) is related to K(p) through 
i^X/u) = pK{p)/{2q — 1). Using the boundary conditions l|13|l. one can determine the 
unknown functions K(p) and C(/i): 

^ r (r7-T + FT^ ri)AiM ) (18) 



C(/x) = - ( J l (m) - J fl (/x) r w - 2(1 - r*)Ai(^) ) (19) 



JV 



9 V 2 « ,-, 

„jV+l\-l 



with B = (1 — r + )~ . Note that the case 5=1, that was temporary excluded from 
the calculation, is recovered by taking the limit q — > 1. A similar calculation can be 
performed in the case q = 1/2 to give: 



K/x) = 2^>' " + jA{n) + 2J4/x) 1 < j < N (20) 



with 

A{n) 



i=i 



N + l 



N 

Mfi) - J L {fi) - ^,{N + 1 - i)Ai(/i) 



(21) 



As mentioned in the above derivation, the form |J5J of the rates <^i(/u|pi) and Ai(p\pi), 
as well as the conditions given in Eqs. I|17|l and l|2UI) are a priori only sufficient 
conditions for the factorization of the distribution P st (pi, Pn)- Yet, let us recall 
that for the general class of mass transport models on a ring geometry studied in 
where dissipation is absent, the form J3J of the transport rate has been shown 
to be a necessary and sufficient condition for factorization, in the case of continuous 
time dynamics. So it may be plausible that this form is also necessary in the present 
model. Let us also note that an explicit test of the form © of the rate functions has 
been proposed in This test should also apply to the present situation so that, 

given two functions <pi(lAPi) an d Aj(//|/5j), one should be able to check whether they 
are of the form J2J), and to determine the corresponding functions 4>i{p), Aj(/x) and 
fi{p)- 
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Site-independent rates 

A situation of physical interest is when the dynamical rules are site-independent 
whereas, due to the injection and dissipation mechanisms, the one-site steady-state 
distribution may depend on the site considered. An example of such a situation has 
been given in |20j. in the fully biased case 5=1. In what follows, we wish to study 
the possibility to find a factorized steady-state distribution with symmetric transport, 
that is for q = 1/2, and with site-independent rates. The transition rates given in 
Eq. © are independent of the site if the functions <j)j (p) , Aj (/x) and fj (p) satisfy 

= /ii(m) e" A ^, AiOu) = MP) e" AjA \ Ij{p) = 9{p) ^ P (22) 
Besides, from Eq. Ijl5|l . the rates have to satisfy 

^•_ 1 (/i) + ^^ +1 (/i)-^(/i)=A J -( / i), l<j<N (23) 

Let us introduce a coordinate x = j/L, so that < x < 1. For the sake of simplicity, 
we shall consider the limit of large system size, and assume that the transition rates 
vary slowly as a function of the site index, in the sense that there exists continuous 
functions (f>(/i,x) and A(p,x) of the variable x, such that 

fa(p) = fo>3L)> Aj-(jj) = -^A(p,jL). (24) 

The 1/L 2 factor has been included so that Eq. (|23l admits a consistent continuous 
limit. From Eq. I|22l) . this means that Xj must be replaced by a continuous function 
\(x). Accordingly, 0(/i,x) and A(p,x) satisfy the following equation 

^(ji,x)=2&(p,x) (25) 
which may be rewritten as 

M A'(*) 2 -A"(,) = ^4 (26) 

As the r.h.s. of this last equation does not depend on x, the only possibility is 
that A(a;) is linear in x, so that the l.h.s. is also a constant. As a result, it turns 
out that there is no factorized solution that would be symmetric in x with respect 
to x = 1/2. This suggests that the present model, with q — 1/2 and symmetric 
injection rates Jl(p) = Jr{(j), does not admit a factorized steady state, at least in 
the continuous limit considered here. The above calculation, with A(x) linear in x, 
rather corresponds to a cascade process, where energy is injected on one boundary 
only, in a way essentially similar to the fully biased case |2(J| . 



A simple physical application 

Coming back to the value q = 1 (fully biased model), one can try to use an 
inhomogeneous version of the model (i.e., with site-dependent rates) in order to 
describe schematically a turbulent cascade -note that only a homogeneous case was 
studied in [20]. To this aim, we assume that the lattice sites n may be thought 
of as successive wavenumbers k n = nki, and that energy is injected on site n = 1 
(corresponding to the largest wavelength) with a rate Jl{p) = Jo e _/3AI , whereas there 
is no energy transfer starting from site n — N, that is Jr(h) = and <ppf(fi\p) = 0. 
Note that the energy flux entering the system, which in steady state equals the total 
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rate of energy dissipation, is given by e = Jq//3 2 . To proceed further, the properties 
of the transition rates have to be specified. For simplicity, we assume /„ (p) to take an 
exponential form, namely f n (p) = e~ anP . A characteristic feature of the turbulence 
phenomena is that the energy dissipation rate fl n on site n is of the form fl n = vk 2 (p n ), 
where v is a viscosity. We thus compute f2„ in the present model, yielding: 



1 



p 



0„ = — - / dpf n (p) / dppA n (p\p) = / dpp,A n (p) (27) 

Jo Jo Jo 

with Z n = J °° dpf n (p), and where the last equality is obtained by permuting the 
integrals and using the form ©. A second constraint arises from Eqs. (|17|l to (|19l) 
encoding the factorization condition for the probability distribution, which reads: 

N 

J L (p) = Y / ^n(p) (28) 

71=1 

A simple form for A n {p\p) that matches the above constraints is 

A n (p) = ^I__./? 2 e -^ (29) 
on condition that the following self-consistency equation is satisfied: 

Then, once A„(/z) is known, one can determine 4* n {p) through Eqs. I|17() to (|19|l . 

In the following, we take the infinite N limit and consider the small viscosity 
regime in which developed turbulence is usually observed. To determine the energy 
spectrum, i.e., the values of a n , we require that the energy flux e does not depend on 
the viscosity in the limit v — > 0. In other words, the a„'s must be such that the sum in 
Eq. (|30|) is independent of v. Assuming that k\ is small and that a n depends smoothly 
on n, one can replace the above sum by an integral over the wavenumber k = kin, 
introducing the energy density E{k) at wavenumber k through E(k n ) = l/(A)ia rl ): 
00 ,,1.2 2 00 r oo 

Y,^^ = ^vklkiE{k n )~ / dkvk 2 E{k) (31) 

n=l n=I 

Making the assumption that the energy p n has the dimension of a velocity square, 
it is well-known from Kolmogorov's K41 theory |21j that a dimensionally consistent 
solution for E(k) that yields a non-vanishing energy flux for small v is given by 

3 ~ 

E(k)=eik-ig(vk), V = (y) * (32) 

where rj is the dissipation scale, and g(x) is a scaling function with a finite limit when 
x — > and such that dxx x ^g{x) = 1. So our model allows one to construct a 
stochastic model for turbulence obeying the K41 phenomenology. Its main drawback 
is that due to the factorization condition, the energy spectrum must be somehow 
included by hand in the stochastic rules, whereas one may rather wish this spectrum 
to arise as a solution of the model, by including in the microscopic rules only some 
basic physical principles. In this respect, it may be of great interest to be able to go 
beyond the factorized steady-state distribution, and one might try to solve this model 
with a more general matrix product ansatz. 
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Conclusion 

In the present note, we have introduced a class of exactly solvable stochastic models, 
which include both transport and dissipation phenomena. The factorization condition 
for the probability distribution, as well as an explicit solution for the latter, have been 
derived. Some specific models within the class may used to illustrate the appearance of 
(generalised) Gumbel distributions in complex systems [201 > or to model, in a schematic 
way, the cascade phenomenon in turbulent flows. More generally, one might hope that 
this class of models could help to understand some fundamental issues concerning the 
statistical behaviour of dissipative systems with a large number of degrees of freedom. 
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